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(T^ ; ABSTRACT 



Quenched chiral perturbation theory is extended to include heavy-light mesons. Non-analytic 
corrections to the decay constants, Isgur-Wise function and masses and mixing of heavy mesons 
are then computed. The results are used to estimate the error due to quenching in lattice 
^1^ , computations of these quantities. For reasonable choices of parameters, it is found that quench- 

Q^ ■ ing has a strong effect on fsj fB, reducing it by as much as 28%. The errors are essentially 

negligible for the Isgur-Wise function and the mixing parameter. 

Qh! (To Appear in Physical Review D) 

^' PACS: 12.38.Gc, 12.39.Fe, 12.39.Hg 
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1 Introduction 

Lattice simulations of hadron properties have made great progress in recent years and there is hope 
that they will soon yield accurate "measurements" of quantities that are difficult or impossible to 
access experimentally, such as the kaon mixing parameter Bk and the B and D decay constants 
/b and fo, play an important role in the phenomenology of the standard model. This progress 
has come not only through improvements in computer speed and algorithms but also through better 
understanding of errors. One systematic error present in most calculations is that arising from the use 
of the quenched (or valence) approximation, in which disconnected fermions loops are neglected. For 
heavy quarks, i.e. those with masses well above the QCD scale, such as the b and c, the decoupling 
theorem ensures that quark loops can be accounted for by suitable adjustments of the coupling 
constants. But for lighter quarks, with masses below the QCD scale, it is expected that quenching 
will change not only the short-distance but also the long-distance properties of the theory; these 
latter changes are much more difficult to quantify. 

A straightforward way to study this error is to perform simulations with dynamical fermions and 
compare the results to similar calculations in the quenched approximation. However, this is still a 
complicated undertaking because calculations with dynamical fermions are performed with larger 
lattice spacings and heavier quarks. This of course increases the errors and makes it more difficult 
to isolate the effect of quenching. For example, the study of fs in Ref. |] sees little effect due to 
quenching, but the interpretation is difficult because of the large lattice spacing used. 

Another approach to understanding the error is to study how quenched QCD differs from full 
QCD in the continuum. That is, one compares the quenched and full QCD predictions for a given 
quantity. Then, to the extent that lattice calculations reproduce the continuum theory, the difference 
between the two predictions gives an indication of the error due to quenching. This analytic approach 
was initiated by Morclg, who studied how chiral logarithms differ in the two theories. It was 
extended by Sharpep|, who developed a diagrammatic analysis, and later Bernard and GoltermanB 
formulated quenched chiral perturbation theory to discuss these logarithms in a systematic way. 
Corrections to light meson decay constants and masses, B^ and recently baryon massesp] have 
been studied using these techniques. 

In this paper I will extend quenched chiral perturbation theory to include heavy-light mesons. 
This enables one to study the effect of quenching on lattice studies of these mesons. The paper is 
organized as follows. In section 2, I review the combination of chiral and heavy quark symmetries. I 
continue the review by showing how chiral perturbation theory can be formulated for the quenched 
approximation to QCD. Finally I show how to extend this to include heavy mesons. In section 3, 

1 compute loop corrections to the heavy meson decay constants and mass splittings, the mixing 
parameter Bb and the Isgur-Wise function ^. In section 4, following a discussion of the parameters 
of the theory, the results are investigated numerically. In section 5 I conclude and comment on 
possibilities for future study. An appendix collects results for the renormalized couplings. 

2 Quenched Chiral Perturbation Theory and the Inclusion 
of Heavy Mesons 

2.1 Chiral Theories 

To lowest order in the chiral expansion, the self-interactions of the light mesons are described by the 
Lagrangian 

C = ^ [tr(a^Sa^I]t) + 2^0 tr(A/S + ME^)] (1) 

8 

with S = ^^ and 

^ = e*'^^^)/^ (2) 



where the Hght mesons are grouped into the usual matrix 

V2^ ^ 76 



i=7rO + 4^7? 7r+ X^ 



-71'^° + Ta'^ ^" . (3) 

The normahzation is such that /^ = 128 MeV. Under SU{2>)l x S'[/(3)fl, ^ transforms as 

^^L£,U^ = U^Rl (4) 

This equation implicitly defines f/ as a function of L, R, and ^. The quark mass matri3i|^ M is given 
by 

(rriu \ 

ma . (5) 

nis J 

For purposes of determining the allowed form of the Lagrangian, M is given the "spurion" transfor- 
mation rule 

M h-^ LMR\ (6) 

so it is convenient to define the quantities 



M+ 



i(CtMCt±^MO, (7) 

which transform as 

M± ^ UM±U'^. (8) 

At leading order in the \/M expansion, strong interactions of B and B* mesons are governed by 
the chiral Lagrangian [|6| 

£ = - tri3 [Ha{v)w ■ DbaHbiv)] + g tm [Ha{v)Hb{v)Aal5] ■ (9) 

The B and B* fields are incorporated into the 4x4 matrix Ha which conveniently encodes the heavy 
quark spin symmetry: 

Ha - i(l+^)[B:%-^a75], (10) 

Ha = 1°HIj\ (U) 

Here v'^ is the four-velocity of the heavy meson, the index a runs over the light quark flavors, u, d, 
s and the subscript "D" indicates that the trace is taken only over Dirac indices. Henceforth I will 
drop explicit reference to the heavy meson velocity. Under SU{3)l x S'f/(3)_R, H transforms as 

H^HU''. (12) 

The light mesons enter the heavy meson Lagrangian Eq. (0) through the quantities: 

V^ = \{^d^S} + ed^^), (13) 



^ There should be no confusion when M is also used to denote a generic heavy meson mass. 



It follows from the definitions that under SU{3)l x SU{3)u 

V^ ^ UV^U^ + iUd^U\ A^ ^ UA^U\ (15) 

while the covariant derivative transforms as 

Df.X^UDf.XUl (16) 

Finally, the left-handed current which mediates the decay B ^ liy is represented by 

J^ = zatrn[r'^Hi,^l] (17) 

where F^ — 7^i ~ 7^(1 — 75)/2- At lowest order the decay constants are related (in my normaliza- 
tion) by fs = a/y/Ms, Jb* = ay' Mb- 

2.2 Quenched QCD 

In the quenched approximation to QCD, the determinant which arises in the functional integral 
when the quark fields are integrated out, is omitted. This can be implemented in a formal way by 
introducing for each quark qa a "ghost" partner (ja with the same mass, but bosonic statistics, so 
that the ghost determinant cancels the quark determinant 0] . The Lagrangian is then 

-Cqucnchod = ^QaiP + ma)qa +^ QaiP + ma)qa- (18) 

a a 

Classically, when the masses vanish, the quenched Fagrangian ( p^ ) is invariant under the graded 
group C/(3|3)l X C/(3|3)_r, but at the quantum level the full symmetry is broken by the anomaly^o 
the semi-direct product^ {SU{3\3)l x 5't/(3|3)_R)(s)t/(l). Elements of the graded symmetry group 
are represented by supermatrices (in block form) 

where A and D are matrices composed of even (commuting) elements and B and C are composed 
of odd (anti-commuting) elements. If we assume that chiral symmetry breaks in the usual way, then 
the dynamics of the remaining 18 Nambu-Goldstone bosons and the 18 Nambu-Goldstone fermions 
can be described by an effective chiral Fagrangian, just as for full QCD |^, ||, |^, |[ |^. The meson 
matrix is extended to a supermatrix 

where x ^ QQi X ^ QQ ^^'^ 'i' ^ 99- Note that x ^'^d x^ si's fermionic fields, while and <j) are 
bosonic. Group invariants are formed using the super trace str and super determinant sdei, defined 
as 

str(C/) = tr(yl) - tr(i:>), (21) 

sdet(C/) = exp(strlog(C/)) =det(A-Bi:)-iC)/det(L>). (22) 

The lowest order Fagrangian would then have the same form as Eq. (|lj) above, with obvious notational 
changes. But because the full symmetry group is broken by the anomaly, extra terms are required 
to describe the dynamics of the anomalous field. In full QCD, this anomalous field is the 77', and 
these extra terms can be neglected because the anomaly pushes the mass of the 77' up beyond the 

■^The broken (7(1) is that which acts as 5 — > e^^^'^q, q — > e~*"^6g. 



chiral scale. However, in the quenched theory, because of the absence of disconnected quark loops, 
this decoupling does not occur: the super-ry' remains in the theory and the extra terms must be 
included. To lowest order, the complete Lagrangian is then 

Cbg - ^ [str(a^sa^st) + 4mo Mm+)\ + ^a^<f>oa^*o - ^*^ (23) 

o Z Z 



with 



$0 = ^str$ = ^(r/'-ry'), (24) 



M - i" t). (25) 



V3 V2 

( M ^ 
I M 

and A^± defined analogously to M±. 

The propagators that are derived from this Lagrangian are the ordinary ones, except for the 
flavor-neutral mesons, for which the non-decoupling of $o leads to a curious double-pole structure. 
For these mesons it is convenient to use a basis C/i, corresponding to uu, dd and so on, including the 
ghost quark counterparts. Then the propagator takes the form 






^V — „2 71^2 "•" /„2 ll^2^^„2 nj2\ y^^l 



where e = (1, 1, 1, —1, —1, —1) and Mf = 2/io'Tii. It is convenient to treat the second term in the 
propagator as a new vertex, the so-called hairpin, with the rule that it can be inserted only once on 
a given meson line. 

Heavy mesons can be incorporated into this framework by adding to H extra fields B and B* 
derived from the heavy fields B and B* by replacing the light quark with a ghost quark. It is 
necessary to include in the Lagrangian vertices which couple $o to H. Symmetry requires that this 
coupling occur through str(j4^), which no longer vanishes. Including also explicit SU{3) breaking 
terms, the Lagrangian is 

C = -tiD [HJv ■ DbaHb] + g tiD [HaHb Aalb] + 7 trp [HaHaJt,!^] str(A'') 

+ 2Ai tri5 [HaHb] {M + )ba + h tn, [Haiv ■ DbcHb] {M + )ca 

+ fe iro '\Haiv ■ DbaHb] str(7W + ). (27) 

The B and B* propagators are i^^ and ^ 2~^^^" , respectively; the ghost mesons have the same 
propagators as their real counterparts. To the same order, the current is given by 

J„^ = mtrz5[r^i/4j+za«itrz5[r^i/4j(M+),6 + m«2trz5[r^i/fc4]str(M+). (28) 

In the sequel, the terms proportional to mq will be loosely referred to as counter-terms because they 
are required to absorb the divergences encountered in loop calculations. In addition, the presence of 
the additional mass scale toq means there will new divergences (not found in the unquenched theory) 
proportional to it. For completeness, the divergent portions of the counter-terms can be found in 
the appendix. 

At this point let mc pause to note a few peculiarities of the theory just formulated. First, while 
the symmetry allows terms involving str(Al-|_), they do not contribute at tree or one- loop level to 
any of the quantities I will consider (although they do contribute in the unquenched theory). Second, 
the loop structure of the quenched theory is rather odd. Because the heavy mesons contain only one 
light quark and there are no disconnected quark loops, none of the meson loops involve any flavor 
changing vertices. Consequently, the loop corrections for a generic heavy meson Bq containing the 
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Figure 1: The diagrams which contribute to the heavy meson self energy. Solid lines represent heavy 
mesons, dashed lines represent light mesons and the cross represents an insertion of the "hairpin" 
vertex. 

light quark q will be a function of Mq alone. The three-flavor theory is then just three copies of a 
single-flavor theory. This tends to heighten the difference between the full and quenched theories 
because not only are the virtual quarks lost, but the "averaging" effect arising from the interaction 
with light mesons of different mass is also lost. This is in contrast with the situation for light 
mesons in quenched QCD, where a kaon has loop corrections involving dd, ss and ds mesons, each 
(potentially) having a different mass. It is the cancelation between these different meson loops, for 
example, which is behind Sharpe's result 0] that Bk is the same in the full and quenched theories 
when TTis = rud- 

3 Loop Corrections 

3.1 Loop Integrals 

There are several loop integrals which will be encountered. Two of these integrals are shown below. 
The first (here Z? = 4 - 2e, 1/e = l/e + log47r - 7e + 1). 

'^"P 1 ~™^1 . 1 _2,_/ 2/„2 



arises from light meson tadpoles, while the heavy-light loops require 



m^ log(TO7/z-^) , (29) 



J^''(to,A) 



(27r)^ (p2-m2)(p.'y- A) 

{(m^ - |A2)i + (|m2 - l^A^) + Mm, A)}cr 
3 e 3 9 



IGtt 

-|-{(2m2 - ^A2)i + ^(7A2 - 3m2) + J2(m, A)}«''«'' 



(30) 



The remaining integrals can be obtained by differentiation with respect to m^, which will be denoted 
with a prime. The definitions of the functions Ji and J2 can be found in Rcf. Bl For my purpose I 
need only the limiting values 

J(m,0) = ^m^ (31) 



dJ{m,0) 
dA 



-m'^log{m^/^l'^), (32) 



where I have defined J{m, A) = AJi(to, A). 

The graphs which contribute to the self-energy arc shown in Fig. m. In diagram la, the ghost 
mesons will cancel the contribution from the real mesons unless one of the vertices involves the 
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Figure 2: Corrections to the weak current vertex. The box represents an insertion of the weak 
current. 



singlet field. Combining this with the contribution of the hairpin vertex diagram lb, I obtain 
6i 



iS(v • k) 



(257 - I g^aQ)J{Md, -vk) + - g^ml - aoMJ)J'{Md, -v ■ k) + . . 



The terms not shown are analytic in Aid and can be obtained from Eq. (pOh above 



(33) 



3.2 Wavefunction Renormalization and Decay Constants 

The wavefunction renormalization constants are obtained by differentiating the self-energy with 
respect to 2u • fc and evaluating on-shell. I find 



1 + 3g^ Hq + 6(57 - - g'^ao) Hd + kiirid. 



(34) 



Here and below, it is convenient to adopt the definitions 



A*P = 



Utt^P "'°' //2 " -'^ 16^2 J2 



"^°/'log(^), 



M 






(35) 



Loop corrections to the left-handed current vertex arise from the diagrams of Fig. g. It is easy to 
see that the diagram Fig. ga vanishes: the loop integral must be proportional to w^, which will 
vanish when contracted with the projection operator in the numerator of the B* propagator. The 
remaining tadpole graph Fig. Ep yields 



iav^i~2aoMi+ml)/3^^^^Ml^_ 



167r2/ 



2 f 2 



M 



(36) 



The final results for the decay constants are then found by combining the wavefunction and vertex 
corrections: 



MbIb = a 



1- 2i^ + 3g^)f^q - (357- (l + 35^)Y)Aid + Kimrf 



(37) 



In contrast, the results in the full theory arepG, pj 



V^^/b. 



Vm^Jb^ = 



1 - -(1 -I- 35^) f - Mtt + MK + T M»7 ) + «lWd + K2{ms + 2md) 
1 - -(1 + 3.9^) f 2 ^K + - ^^ j + Kinis + K2{ms + 2md) 



(38) 
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Figure 3: The corrections to B — B mixing. The double box represents an insertion of the mixing 
operator. Corresponding diagrams with hairpin vertices are not shown. 

3.3 Masses 

The correction to the mass is obtained by evaluating the self-energy on-shell and removing the 
wavefunction renormalization constant (though to the order I am working it does not contribute). 
Defining 

Mb =Mb-^A + SMb (39) 

where Mb is the spin- averaged mass in the chiral limit, A is the hyperfine splitting and SMb is the 
light-quark dependent contribution to the mass, I find 

6Mbo = 2Ai m, - ^^ i^g' ^M, + (2^7 - 5<?^|^)Afj) . (40) 

while in the unqucnched theory pj|, O] 

2 

SMb^ - SMb, - 2Ai(m, - nid) - T^^{~Sml + 2m\ + ml). (4f) 

3.4 Mixing 

The constant Bb^ is defined as the ratio 

^ _ (BairLl^bLqlrbLlB,) _ {Ba\rLl^bLqlrbL\B,) 



l{Ba\qllMO){0\rLrbL\Ba) |/K„ Bb. 



(42) 



As shown by Grinstein and collaborators [l0| , in the effective theory the operator 

TLlf^bLqirbL (43) 

is represented by 



4/3 trn 



i^H^'^rj^L] tYD Im^'^rrL , (44) 



which is essentially just the square of the left-handed current. The one-loop corrections to this 
operator are shown in Fig. (0). There are two types of tadpoles which arise from the operator 
Eq. (E|): those where each ^ is expanded to O(^) (Fig. ||a) and those where only one of the ^'s is 



expanded (Fig. 0b). The latter tadpoles just renormalize fs and will cancel in the ratio for Bb- 
Thus it is only necessary to consider the former. I find 



Bb, = 4/3 



1 - (1 - 35 ) ( (1 - -ao) /"d + A*9 ) + 637 fJ-d + Pirrid 



(45) 



which should be compared with the unquenched results |lO[ {(3i and P2 are additional counter-terms) 

2 



Bb = 4/3 
Bb. = 4/3 



1 - (1 - 35 )- A*») + f^imd + Piirris + 2md) 

1 - (1 - 35^) f ^ Mtt + ^ A'r, ) + Pims + /32(m-s + 2TO<i) 



(46) 
(47) 



The reader will note that in contrast to the earlier results, Bb has true chiral logarithms even in 
the absence of the singlet coupling 7 and the kinetic coupling ao- A similar phenomenon occurs in 
Bk, as shown by SharpeO]. The reason is that flavor conservation allows disconnected quark loops 
to appear only in the guise of the 77', so that even in full QCD they do not contribute. 

3.5 Isgur-Wise Function 

The heavy quark current which mediates the decay B — > Dev is represented at leading order by 

U^)iro[H''-\v')l,LH^^\v)]. (48) 

where S,o{'w) (w = v ■ v') is the leading-order Isgur-Wise function. In the full theory, the leading 
corrections are[|ll[, |l3| (here the counter-terms 771 and 772 are functions of w) 



SU.Q 



Co(^) 



3 1 

1 + 2g^{r{w) - 1) ( - A^TT + A*K + g M>) 



-I- 2(r(w) - 1) rii{w)md + ri2{w){ms + 27nd) 



?o(w) 



l + 2g\r{w)-l) (2AiA' + -M„ 



where 



The quenched results take the by-now-expected form 



-I- 2(r(w) - 1) f rii{w)ms -f r]2{w){ms + 2md) 
r{'w) = - log(w + \/w'^ — 1). 



U{w) =Co(w) 



l + 2(r(w;) - 1) [g^ i-J.q + (2.g7 - -g'^ao) ^d + Viiw)md 



(49) 
(50) 

(51) 



4 Discussion and Numeric Results 

To obtain numeric values it is necessary to know the values of the various couplings which enter 
the Lagrangian. Combining data on the D* width and branching fractions M, 15|, Amundson et al. 
Rcf. dig ) obtained the constraint 0.1 < g^ < 0.5. The spread is caused by the uncertainty in the 
branching fraction BR{D*^ -^ D^l)\ taking the central value yields g ~ 0.5. QCD sum rules fT^ 
and relativistic quark models[18| favor a smaller value, g ^ 1/3. Given this uncertainty, I will show 



Table 1: Quenched and full heavy-light quantities. Only the non-analytic part of the mass difference 
is shown. 



Quantity 


Unqucnchcd 


mo = 750 MeV 


mo = 1100 MeV 


fsJfB,-! 


0.074(1 + 3.92) 


-0.11(1 + 352) 


-0.23(1 + 352) 


BbJBb,-1 


0.052(1 - 3.g2) 


-0.11(1-352) 


-0.36(1-352) 


^s/^d - 1 


0.05952 


-0.086 52 


-0.18 52 


SMb, - SMBa 


-450g2MeV 


-340 52 MeV 


-740 52 MeV 



results for different values of the coupling. There is no information on the coupling 7, but 1/Nc 
arguments suggest that it is small.0 They also suggest that ao is small; direct evidence from 7] — rj' 
mixing confirms this. Consequently, I will take both 7 and ao to vanishjj The maximum value of 
w in the decay B -^ Dlv is about 1.8, so I will use r(1.8) = 0.76 when evaluating ^. Finally, I will 
choose ^JL—l GeV. 

There are several ways to determine vtiq, each giving a different result. The Witten-Veneziano 
large Nc formula] 22l rri^ = rri^, + rri^ — 2m\ gives mo « 852 MeV, while from the r] — rj' mass 
sphtting Sharpcja] estimated mo ~ 900 MeV. It has also been computed directly on the lattice. 
Early attempts pi | found vtiq « 570 — 920 MeV, but with limited statistics and a strong dependence 
on the lattice spacing. Recently, a more accurate computation has been performed. Kuramashi 
et al.^^ extracted mo by comparing the one and two loop contributions to the 77' propagator; 
they found mo — 751(39) MeV. Using the U{1) Ward identity relation mg = &xl fii with x the 
topological susceptibihty, the same group found mo = 1146(67) with x and /^ obtained on the same 
lattice. They attributed this larger result to contamination from extra terms in the Ward identity 
induced by the use of Wilson fermions. I will choose m-o = 750 MeV, but will also show some results 
for mo = 1100 MeV. 

For an honest calculation, it is also necessary to specify the 0{mq) counter-terms. But it is an 
unfortunate fact that there is little to constrain them, save the general expectation that their natural 
scale is A^ « 47r/. A common practice when confronted with this situation to assume that the 
counter-terms are overshadowed by the logarithmic contributions. Another approach is to reduce the 
dependence on these unknown terms by taking appropriate ratios. Fortunately, since it is expected 
that the coefficients in the chiral expansion should be (almost) the same in the quenched and full 
theories, some of the counter-terms will cancel when the predictions of the two theories are subtracted 
to compute the error. In particular, the errors in /b^/Zb^, Cs/Cdi Bb^/Bb^ and SMb^ — SMb^ are 
free of counter-terms. Results for these quantities are shown in Tabic |l|; in order to illustrate the 
m,o dependence, the quenched results are shown at both m-o ~ 750 MeV and mo = 1100 MeV. The 
ratios in the quenched theory are computed by substituting Ms = ^2m.\^ — m^ = 680 MeV and 
Md = m-K. Concentrating on the results at rriQ — 750 MeV, one sees that the corrections to the 
ratios are similar in magnitude but opposite in sign to those in the full theory. This is a result of 
the fact that the quenched logarithms diverge in the chiral limit. Notice that the corrections for the 
mass splittings threaten to be larger than the splittings themselves unless 5 is small. This suggests 
either a large cancelation occurs with the leading Ai term or higher-order corrections are important. 
Either solution casts doubt on the reliability of the error estimates in this case. 

While the results in Table |l| suggest large quenching errors — particularly in JbJ Ib — it is 
likely that the error in actual simulations will be less. The reason is the following. Currently, most 

^ The same argument implies a suppression of the singlet coupling to the nucleon. This is confirmed in a phe- 
nomenological study by Hatsuda|19[, who found g^iNN ^ t-1) which should be compared with Qt^mn = 13.4. 
^ Even if 7 is as large as g/3, I find it changes the quenched results by only 5% or so. 
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Figure 4: Quenched and unquenched corrections to /s as a function of meson mass. Also shown is 
a hnear (in the quark, not the meson mass) extrapolation of the quenched fs- The horizontal lines 
are the predictions for fs^ and fs in the full theory. 

simulations are performed with quark masses corresponding to pion masses in the range 400 ^ M-^ < 
1000 MeV. The results are then extrapolated linearly (in the quark mass) to the chiral limit and 
the physical ttItt ■ Due to the familiar property of the logarithm, quenched loop corrections change 
as much in the interval 140 < Md < 350 MeV as they do in the interval 350 < Md < 1000 MeV. 
Consequently, in the mass range covered by lattice simulations the quenched logarithm appears 
linear. This can be seen in Fig. Q, where both the "true" quenched and linearly extrapolated 
predictions for /b are shown. Clearly, the two cannot be distinguished for masses greater than 
say 300 MeV, but the extrapolated result underestimates the "true" behavior by more than 10% at 
Md = rriT^. While the primary motivation for this extrapolation is the desire to efficiently invert the 
quark propagators, it has the side-effect of reducing the quenching errors. Moreover, it is the correct 
thing to do, since the goal is to describe unquenched QCD and quenched ChPT clearly fails to do 
this in the chiral limit. Table || compares the two methods of computing the error at different values 
of the coupling. One sees that the error is substantially reduced by the extrapolation. The errors 
in Bb and ^ are relatively small to start with and become negligible when extrapolated. However, 
even with the extrapolation the error in fs is larger then one might have hoped. Note also that 
/Ss //b is smaller in the quenched theory. 

In fact, the size of the error in /b is easy to understand. For the Isgur-Wise function, the error 
is small because the corrections themselves are small. Conversely, the corrections to /b are large 
and so the error is large. Moreover, they are driven by the tadpole terms, which remain large even 
if the coupling g vanishes. The tadpoles, however, do not depend on the heavy quark mass, so it 
should be possible to eliminate them by studying the 1/M corrections p3[ . 

Some additional understanding of the differences between the quenched and dynamical theories 
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Table 2: The errors for different choices of the coupling g. Results are shown using both the exact 
and extrapolated quenched predictions. An error is negative when the quenched quantity is smaller 
than the full. 



Errors 


5 = 0.7 


Exact 
0.5 


0.33 


Extrapolated 
0.7 0.5 


0.33 


IbJIb, 


-0.45 


-0.32 


-0.24 


-0.28 


-0.20 


-0.15 


BbJBb, 


0.078 


-0.041 


-0.11 


0.012 


-0.0062 - 


-0.017 


^s/U 


-0.072 


-0.036 


-0.016 


-0.045 


-0.023 - 


-0.010 


Mb, - Mb^ 


49MeV 


25 MeV 


11 MeV 


- 


- 


- 



may be gained by comparing them in the mass range probed on the lattice. For this it may be better 
to consider a two quark theory with degenerate masses (rather than the full three- flavor theory), 
since it is closer to the type of theory studied in unquenched simulations. In Fig. ^ the predictions 
for Ib E^re shown. Here I have neglected the unknown counter-terms, though it is clear from the 
graph that there must be a positive term of 0{md) since lattice simulations find that Jb increases 
with the light quark mass. There are two general features that should be noted in Fig. ^. First, 
that in both the quenched and two-flavor theories, /s is less than both fs^ and Jb,- This may be 
attributed to the fact that the full theory has more mesons contributing to loop corrections. The 
second observation is that the gap between the two-flavor and quenched corrections grows as M^ 
increases toward the point Mg = -y/2m^ — rn^. Thus, the fact that quenching decreases the ratio 
IbJIb^ is due to the different nature of the quenched logarithm. Finally, Table ^ compares the 
quenched and two-flavor predictions at a representative mass of Md = 600 MeV (again neglecting 
counter-terms). It can be seen that the errors are comparable to those found in the extrapolated 



Table 3: Quenched and unquenched results at Md {— TOtt) — 600 MeV. 



Quantity 


TOO = 750 MeV 


Unquenched {Nf = 2) 


VMfB/a - 1 


0.031(1 + 3.g2) 


0.10(1 + 3^2) 


Sb/4/3 - 1 


0.20(1 - 3g^) 


0.069(1 -3^2) 


Uio - 1 


0.025 g2 


0.083 ,g2 



ratios. 



5 Conclusions 

1 have included heavy mesons into the framework of quenched chiral perturbation theory and used it 
to study the error arising from the use of the quenched approximation in lattice studies of heavy-light 
mesons. These lattice studies are important for the phenomenology of the standard model. Because 
estimates of the error depend on the as-yet unmeasured value of the the B*BTr coupling g, results 
were shown for several values of g in the allowed range. It was seen that the errors in Bb and S.{w) 
were negligible. However, the error in /b was surprisingly large, more than 15% in the best case. It 
was observed that the large error follows from the large corrections present in both theories. These 
large corrections were traced to the tadpole corrections, and it was suggested that they might be 
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eliminated by studying the \/M corrections to the theory. Indeed, Grinstein[E3, M has advocated 

f / f — 

doing just that in the continuum by studying the double ratio f ° / f ■ This is done for the quenched 

theory in a forthcoming work ||23| . 

A general conclusion that can be drawn from the quenched chiral calculation is that quenching 
tends to decrease the ratio JdsI Id- This is in agreement with the one unquenched simulation pi, 
which found foslfo = 1-34, a value larger than that typically found in quenched calculations. This 
may have implications for reconciling lattice predictions of f^^ with the recent CLEO measurement P5| 

In the future it would be interesting to study heavy baryons containing two light quarks within 
this framework. Because of the presence of two light quarks, the quenched theory will be less trivial 
and the loop corrections will have a more complicated flavor structure, more closely resembling that 
of the light mesons. It would also be useful to to move beyond 1/Nc arguments for the magnitude 
of 7. It appears that it could be calculated within QCD sum rules using the same techniques that 
have recently been apphed to g]17|. 
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A Renormalization Constants 

Within the context of dimensional rcgularization, the singularities of the effective Lagrangian are 
customarily described in terms of the parameter L(/i) which contains the singularity at D = 4 (recall 
D = 4-2e, l/e= 1/e + log47r - 7e + 1): 



L{^x) 



1 



IGtt- 



rM 



(A.l) 



It is convenient to decompose an arbitrary coupling k as k = fc'^(/i) + k L{^). 

To render the Lagrangian Eq. ( |27| ) finite, it is necessary to add the counter-term 



In addition, ki must be taken to be 



P 



ki =Q{g^- -g'^ao) 



2 ^i^/fo 
"" P' 



The current Eq. (E8) is renormalized with 



3 



-1 = -((1 + 35^)^ -3.97) 9^. 



2^0 
P 



and in addition a must be rescaled: 

a' (/i) = a 



1 + ^(1 + 3<?^)^L(m) 



The description oi B — B mixing requires the counter-term (no sum on a) 

iPf3ltlD [(CS^^'VlM^l trz5 [(ei/(^')"7^^1 {M + )aa, 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 
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and the couplings must be 



p^Ui) = /3(l + (l-3.g2)!^i0,)^, (A.7) 



/3i = ((l-3ff2)(l--ao)-657)^. (A.8 



Finally, the Lagrangian for b ^ c transitions needs the countcr-tcrni 

r,iw){riw) - 1) tTD[H['\v'h^LHi''\v)]{M + )ba, (A.9) 



with 

1 2 ^ 2Aio 

and the rescaled coupling 



rjiw)^-4{gj--g^ao)^ (A.IO) 



^S(i«, m) - CoH ( 1 + 2(rH - 1) .g^^L(/i) ) . (A.ll) 
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